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Abstract 


A set  of  constitutive  equations  has  been  developed  to 
represent  nonlinear  anelastlc  behavior,  l.e.  energy  losses  amd 
rate  dependent  moduli  for  geometrically  large,  reversible 
deformations.  These  equations  are  in  incremental  form  and  can 
be  used  to  solve  boundary  value  problems  by  numerical  methods. 
Two  problems  were  considered:  thick  walled  spheres  subjected  to 
cycles  of  internal  pressure,  and  a long  rod  subjected  to  indeed 
axial  velocity  (the  simple  tension  test) . Experiments  were 
carried  out  for  both  conditions  on  specimens  made  of  a synthetic 
rubber  at  three  different  temperatures  (each  under  isothermal 
conditions)  and  at  different  loading  rates.  Response  curves 
were  calculated  based  on  the  constitutive  equations  and  chosen 
values  of  the  material  constants.  Generally  good  agreement  was 

obtained  between  the  predicted  and  experimental  results.  
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Relatively  little  work  has  been  done,  however,  on  formulating 
emelastic  behavior  in  the  nonlinear  range.  This  paper  describes 
I , a procedure  for  representing  nonlinear  anelasticity  which  is 

[ similar  to  one  developed  for  inelastic  (non-reversible)  deforma- 

^ tions  by  Partom  amd  Bodner  [6,7,8].  The  constitutive  equations 

consider  the  increments  of  both  the  reversible  emd  the  non-rever- 
sible deformation  rate  conponents  to  be  functions  of  state  queui- 
tities  and  the  current  geometry.  The  chemges  in  the  state  quan- 
tities with  deformation  are  given  by  suitable  evolutionary 
equations.  In  this  manner,  the  material  response  to  arbitrary 

I 

I 

loading  or  straining  historijes  is  calculeJsle  on  a step  by  step 
basis.  Anelastic  behavior  cam  be  considered  in  a similar  manner  by 
l'  an  anelastic  stress  term,  not  a state  variable,  which  is  a non- 

I linear  function  of  the  total  deformation  rate  and  state  quauitities. 

i Only  reversible  deformations  are  considered  in  the  present  paper 

but  generally  both  smelasticity  and  inelasticity  can  be  treated 
simultaneously. 

Since  the  deformations  are  reversible  and  the  total  stress 
is  considered  to  be  decomposedjle  into  elastic  and  emelastic  com- 
ponents, the  present  formulation  is  essentially  equivalent  to  a 
nonlinear  generalisation  of  the  Kelvln-Voigt  viscoelastic  model. 
Nonlinearity  is  Introduced  in  the  elastic  element  through  the 
use  of  a general  nonlinear  stress-strain  relation  whose  coeffi- 
cients are  derived  from  a strain  energy  function  and  by  the  in- 
clusion of  a finite  strain  measure.  The  anelastic  stress  is  a 
nonlinear  function  of  the  deformation  rate  as  well  as  a function 
of  the  elastic  stress.  One  of  the  consequences  of  this  approach 
i is  that  the  effective  "dynamic  modulus"  and  "damping"  become  un- 

coupled so  their  frequency  dependence  could  be  specified  separ- 
ately idiioh  is  essentially  different  than  linear  viscoelastic 
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theory . 

Lockett  in  [1,2]  classified  nonlinear  viscoelasticity  as 
a consequence  of  various  effects:  "magnitude,"  i.e.  stress  emd 
strain  amplitude,  "interaction,"  i.e.  different  times  of  load 
application,  and  "intermode,"  i.e.  chsmges  in  direction  of  load 
application.  In  his  terminology,  the  present  formulation  atten?)ts 
to  treat  both  magnitude  and  interaction  nonlinearity.  Intermode 
nonlinearity  is  also  inplied  but  only  investigated  insofar  as 
two  different  boundary  value  problems  are  treated.  There  are 
a large  number  of  interaction  effects  which  could  be  examined, 
such  as  superimpoaed  loading  during  creep  or  recovery,  euid 
Lockett  has  developed  a modified  Kelvin-Voigt  model  to  treat 
such  problems  [91.  The  interaction  exeunples  in  this  paper  are 
limited  to  the  frequency  dependence  of  the  load-elongation 
relation  during  loading  and  unloading  cycles. 

A brief  review  of  the  theoretical  development  of  the 
nonlinear  anelastic  formulation  was  presented  in  [10]  where 
some  results  were  also  given  for  thick  walled  hollow  spheres 
subjected  to  loading  and  unloading  cycles  of  internal  pressure 
under  quasi-static  emd  isothermal  conditions.  This  paper  gives 
a more  complete  presentation  of  the  basic  formulation  and 
describes  and 'discusses  the  sphere  results  in  greater  detail. 
These  spheres  were  prepared  from  a synthetic  rubber  and  were 
tested  at  various  frequencies  and  temperatures.  In  addition, 
an  extensive  study  was  made  on  the  same  material  for  the  uniaxial 
loading  condlti<Hi  based  on  the  same  constitutive  equations  and 
material  constants.  Both  analytical  and  experimental  results 
for  this  case  were  obtained  and  are  discussed  in  this  paper. 
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Nonlinear  Anelastic  Model 

A general  set  of  constitutive  equations  to  represent  in- 
elastic and  anelastic  material  behavior  at  finite  strains  was 
initially  proposed  by  Partom  111] . These  were  subsequently 
developed  by  Partom  and  Bodner  in  a series  of  papers  [6,7,8]  to 
apply  to  the  elastic-plastic  euid  viscoplastic  response  of  both 
non-hardening  amd  strain  hardening  materials.  In  this  general 
formulation,  the  deformation  rate. 


(1) 


where  v^  is  the  particle  velocity,  is  assumed  decomposad>le  into 
geometrically  reversible  and  non-reversible  consonants,  d^^  amd 
d^^.  The  stress  is  taken  to  have  an  anelastic  component  0*^  to 
account  for  energy  losses  (hysteresis)  in  reversible  motions  due 
to  internal  non-homogeneities.  The  elastic  stress  is  a state 
quantity.  Under  these  assumptions,  the  rate  of  work  per  unit  of 
current  volume  can  be  expressed  as 

W - (ajj  + ojj) (d®^  + dPj)  (2) 

The  deformation  rate  components  and  the  anelastic  stress 
are  considered  to  be  functlcms  of  state  variables,  e.g.  j > 
and  the  velocity  field.  All  history  and  memory  effects  are 
therefore  included  in  the  current  state  quantities  and  the  geo- 
metry so  that  problems  could  be  solved  on  a direct  incremental 
basis  in  tlSMi.  This  makes  the  formulation  readily  adoptable  to 
the  computer  solution  of  boundary  value  problems. 
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In  the  present  treatment  of  anelastlc  materials,  the  non- 
reversible  deformation  rate  d?^  is  taken  to  be  identically  zero 
although,  in  principle,  both  anelastlc  and  inelastic  effects 
could  be  considered  simultaneously.  The  problem  therefore 
reduces  to  determining  suitable  expressions  for  d^^  (subsequently 
identified  as  dj^j)  and  for  . 

Firstly,  suitedile  constitutive  equations  are  obtained  for 
d^j  as  a function  of  the  elastic  stress  on  the  basis  of  non- 
linear (large  deformation)  elasticity  theory.  Using  the  Almansi 
finite  strain  measure,  the  relation  between  the  strain  rate  and 
the  deformation  rate  is  given  by  [6,  12] 

"ij  - ■*ij  - <3> 

where  denotes  the  total  strain.  A general  expression  for  the 
elastic  "true"  stress-strain  relationship  is,  from  [6]  and  [14] , 

ajj  - Poc5^j  + 

where  the  coefficients  are  obtained  from  strain  energy  functions 
(as  described  in  [<S]).  Eq.  (4)  together  with  the  strain  energy 
functions  and  the  strain  rate  expression,  Eq.  (3) , would  then 
define  the  elastic  constitutive  equations.  An  uncoupled  form  of 
the  strain  energy  function  was  developed  in  [13]  and  used  sub- 
sequently In  [6]  and  [7].  Since  the  details  of  the  development 
and  the  resulting  constitutive  equations  are  given  essentially  in 
[6]  and  [7],  they  will  not  be  repeated  here. 
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Modifications  to  the  previous  work  were  required,  however, 
to  obtain  adequate  representation  in  certain  cases.  In  [6]  and 
[7],  the  expression  for  the  distortional  energy,  which  is  taken 
to  be  uncoupled  from  the  dilatational  energy,  is  written  in  terms 
of  the  Almansi  strain  invariemts,  Ei,  £2,  as 


W = 6yo(E*  - 3E2)/(3  - 2Ei)* 


where  uo  is  the  classical  elastic  shear  modulus  (see  also  [14]). 

Eq.  (5)  is  a particular  case  of  a more  general  form  for  W based 

s 

on  the  shear  energy  being  homogeneous  of  degree  zero  in  the 
principal  extensions,  neunely 


w » r ( ai(3-2Ei)^  + 02 (3-4Eit4E2)  1 
s ^ ^ Pi(3-2Ei)i  + 62 (3-4Ei+4E2)  ^ 


where  C and  the  o' s zmd  g's  are  constants.  In  order  for  the 
shear  energy  to  vanish  for  pure  dilatation,  3oi  + 02  * and 
to  avoid  singularities,  33 1 -i-  62  7^  0.  In  addition,  for  small 
strains  W should  reduce  to  the  classical  value  (-2yoE2).  For 

S * 

ai  set  equal  to  unity  as  a common  factor,  these  conditions  lead 

to  02  “ - 3 euid  32  “ (l-93i)/3  so  that  W reduces  to 

s 


2yo  { Y ~ 


El  - 3Es 

rrp-Ti 


Eq.  (5)  is  obtained  from  Eq.  (7)  by  setting  3a  * 0 (i.e.  3i«l/9) 
which  was  found  adequate  for  strains  up  to  about  IS  to  30%  (de- 
pending on  teat  temperature) . For  strains  as  high  as  50  to  100% , 
it  was  found  that  nonzero  values  of  3a  led  to  better  overall  rep- 
resentaticxi  of  the  experimental  results. 


Another  modification  to  the  strain  energy  function  was  to 
introduce  wecik  coupling  between  the  shear  and  dilatational 
energies.  This  can  be  done  by  considering  the  shear  modulus  y 
to  depend  on  the  density  change  due  to  compressibility,  i.e. 


y = yo (So) ^ 


(8) 


where 


2 

So 


(p/Po)  ^ 


l-2Ei  + 4E2  - 8E3 


(9) 


zmd  (p/po)  is  the  ratio  of  current  to  initial  density  and  the 
E's  are  the  Almansi  strain  invariants.  For  the  material  auid 
conditions  examined,  y could  be  set  equal  to  zero  in  most  cases 
but  not  all.  For  con^leteness,  it  is  noted  that  the  expression 
used  for  the  dilatational  energy,  as  in  [6]  and  [7] , is 


W 


V 


9 Ko 

7 STF=5T 


la 

[(1-So) * 


- ^(1-So) 


(10) 


where  Ko  is  the  bulJc  modulus  and  a and  3 are  constants  that 
determine  the  variation  of  compressibility  with  pressure.  As 
in  the  previous  work,  we  used  a ^ 1 and  3 2. 

The  elastic  shear  modulus  yo  was  teOcen  to  be  a constamt 
Independent  of  temperature  over  the  range  of  tenperatures  in 
which  the  material  behaves  in  a viscoelastic  mamner.  That  is, 
the  shear  modulus  value  in  the  rubbery  region  was  taken  to  be 
the  const2mt  value.  Temperature  variations  were  Introduced  in 
the  parameters  that  make  up  the  time  dependence,  i.e.  in  the  « 
expression  for  the  anelastic  stress,  which  leads  to  an  inter- 
relationship between  temperature  and  time  which  is  generally 
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observed  for  polymeric  materials.  It  was  necessary,  however,  to 
consider  Ko  to  have  some  temperature  dependence  to  correlate  with 
the  experimental  data. 

In  the  present  analysis,  the  anelastic  stress  term  is 
specifically  introduced  to  account  for  energy  losses  under  geo- 
metrically reversible  deformations.  In  addition,  it  should  repre- 
sent other  effects  such  as  rate  dependent  dynamic  moduli  and  traui- 
sient  effects  during  creep  and  stress  relaxation.  Anelasticity  is 
associated  essentially  with  shear  deformations  so  that  the  terms 
in  the  governing  law  should  be  the  deviatoric  components  of  stress 
amd  deformation  rate. 

General  observations  of  anelastic  effects  in  many  materials 
indicate  that  the  dyneunic  modulus  increases  with  increasing 
frequency  but  that  the  damping  factor,  tan6  or  Q~^,  changes  only 
slightly  with  frequency  over  a wide  range  except  near  a relaxa- 
tion point.  The  deunping  factor  is  defined  in  the  usual  meuiner 

tan6  = W^/2iTWg  (11) 

where  amd  are  respectively  the  dissipated  and  the  maximum 
elastic  energy  in  a cycle  per  unit  original  volume.  Because  of 
the  eUsove  mentioned  characteristics  of  real  materials,  the  standard 
viscoelastic  models,  which  specify  an  interrelation  between  the 
frequency  dependence  of  the  dynaunic  modulus  and  the  dancing,  can 
only  be  fitted  to  experimental  data  in  a narrow  frequency  band. 

The  present  formulation  enables  the  damping  and  modulus  changes 
with  frequency  to  be  expressed  Independently  although  some  inter- 
action results  from  the  nonlinearity. 
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Making  the  anelastic  stress  term  independent  of  time  would 
lead  to  frequency  independent  energy  losses.  A proposed  expres- 
sion for  the  energy  loss  term  is 


(12) 


where  f is  a constant/  the  bar  symbol  indicates  the  deviatoric 
con?)onent,  and  is  the  absolute  value  of  the  first  invariant 

of  the  deformation  rate  (rate  of  volume  change) . Eq.  (12)  is 
similar  to  a Coulomb  type  dry  friction  and  implies  that  the 
friction  magnitude  depends  on  the  ratio  of  the  deviatoric  to  the 
dilatational  deformation  rates.  It  is  noted  that  dilatational 
deformations  would  always  be  present  due  to  nonlinearity  emd 
possible  coupling  factors  euid  experiments  have  shown  that  volume 
changes  do  influence  the  energy  loss  factor  (15,161. 

Rate  dependent  material  response  could  be  introduced  by 
making  depend  directly  on  and  the  deformation  rate. 

A proposed  equation  for  this  effect  is 


( 


■ 2ij  ^ exp(-a^/(-D.)‘/'] 

^ n=i 


(13) 


where  g and  a are  constants  auid  Da  is  the  second  invariant  of 
the  deviatoric  component  of  the  deformation  rate.  Since  Da 
is  a negative  quantity,  the  minus  sign  is  Introduced  to  make 
the  term  (-D2)  positive.  The  series  form  of  Eq.  (13)  is 
zuialogous  to  that  used  to  express  linear  viscoelastic  func- 
tions in  order  to  approximate  experimental  results  (Prony 
Series) . Eq.  (13)  corresponds  to  an  additional  increment  to 
the  elastic  stress  whose  magnitude  increases  with 
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deformation  rate.  Using  only  Eq.  (13)  in  conjunction  with  the 
previous  elastic  equations  would  lead  to  a rate  dependent  non- 
linear elasticity  theory,  i.e.  the  material  response  would  be 
rate  dependent  but  fully  reversible  without  energy  losses.  The 
amelastic  effects  of  importance  in  the  present  investigation  are 
therefore  given  by  Eq.  (12) , which  provides  for  rate  independent 
energy  losses,  and  Eq.  (13)  for  frequency  dependent  response. 
Coupling  of  these  effects  will  occur,  however,  because  of  the 
essential  nonlinearity  of  the  formulation.  The  total  expression 
for  the  amelastic  stress  term  for  the  present  problem  is  then 


2*j 


N 


n' 


(14) 


which  is  used  in  conjunction  with  the  elastic  stress -strain  re- 
lations to  give  the  set  of  constitutive  equations  for  anelastlc 
materials  in  the  solution  of  boundary  value  problems.  In  these 
problems,  the  total  stress  in  the  equilibrium  equations  is  the 
sum  of  the  elastic  and  anelastlc  components. 

The  constitutive  aquations  consisting  of  the  elastic  and 
anelastlc  terms  were  applied  to  coirpute  the  quasi-static  response 
of  thick  walled  spherical  shells  made  from  a synthetic  rubber 
subjected  to  cycles  of  internal  pressure  at  various  frequencies 
and  at  three  constant  temperatures.  A description  of  these 
ejqperimsnts  on  spheres  is  given  in  the  next  section  and  a short 
review  of  the  investigation  has  been  presented  [10].  In  that 
application,  the  equilibrium  equation  was  solved  at  each  step 
in  time  at  a number  of  points  along  the  radial  coordinate  ac- 
cording to  the  required  accuracy.  The  state  (stresses,  deforma- 
tions, deformation  rates)  at  the  end  of  each  time  stop  were  used 
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as  the  basis  for  the  calculation  of  the  new  quemtities  at  the 
next  step. 

For  the  uniaxial  straining  case,  the  S2une  constitutive 
equations  were  applied  to  calculate  the  stress  resulting  from 
extension  of  a thin  long  rod  of  the  saune  synthetic  rubber. 

Tests  auid  corresponding  calculations  were  carried  out  for  three 
inposed  axial  velocities  at  the  same  three  temperatures.  The 
calculations  were  initially  based  on  the  same  equations  and 
material  constemts  used  for  the  sphere  problem.  It  was  found, 
however,  that  the  use  of  Eq.  (5)  for  the  shear  strain  energy 
limited  the  predictive  capability  of  the  calculations  for  the 
tensile  straining  problem  to  a maximum  strain  of  ad)out  30%. 

Eq.  (5)  is  derived  from  the  more  general  expression  for  W , 

S 

Eqs.  (6)  and  (7),  by  setting  the  constemt  3a’‘0.  Considering 

amd  using  Eq.  (7)  for  W in  which  (l-9$i)  *38a  led  to  very 

s 

good  fitting  of  the  uniaxial  tests  results  for  strains  up  to 
100%.  This  change  had  very  little  effect  on  the  calculations  for 
the  sphere  problem  over  the  range  of  Interest.  An  auxillairy  pro- 
gram was  used  to  determine  82  based  on  matching  a given  point  on 
the  stress-deformation  curve  f cr  a particular  crosshead  velocity. 

The  Almansi  strains  for  the  uniaxial  case  are  simply 
^ij  “ (1-XJ*)  for  i-j  and  ■ 0 for  ij<j  where  the  are 

the  extension  ratios.  It  is  noted  that  the  terms  Idj^j^l  end  Df 
in  Eq.  (14)  could  be  expressed  as  functions  of  Poisson's  Ratio 
since  dj^j^  ■ end  Da  ■ ^jXa  + Xa^i  - (l/3)dj^j^. 

For  uniaxial  straining,  ^a  ■■  ~ vli,  so  that  d^^j^  - ij(I-2v) 

and  Da  - - (1/3)  d+v)  * . 

The  steps  in  the  numerical  procedure  and  the  associated  bloc)c 


diagram  ere  given  as  follows} 


I 


4.  Using  the  values  of  0i  obtained  in  step  3,  the  axial  stress, 

01,  is  computed  for  a series  of  deformations  during 

loading  and  unloading  using  a double  phase  Newton -Raphs on 
iteration  procedure  and  applying  the  conditions  (a)  and  (c) 
described  above.  The  convergence  is  speeded  up  by  the  use 
of  relaxation  factors  and  using  initial  values  based  on 
results  of  the  preceding  steps. 

5.  Using  the  same  value  of  3i,  step  4 is  repeated  for  other 
Imposed  velocities. 

The  block  diagram  of  the  calculation  procedure  for  uniaxial 
extension  is  given  on  the  following  page. 


I 


! 


i 
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Experimental  Procedure  and  Results 

For  the  thick  walled  sphere  problem,  two  sets  of  spheres 
having  inner-outer  diameters  of  60  - 120  mm  and  75  - 100  mm 
were  fed>ricated  of  a synthetic  rubber.  To  measure  the  inner 
radial  displacement,  a specially  designed  extensometer  was  em- 
bedded in  the  inner  cavity  prior  to  casting  the  sphere  as  a 
single  piece.  Another  extensometer  was  attached  to  the  outer 
surface.  The  spheres  were  subjected  to  rectified  sinusoidal 
cycles  of  positive  internal  pressure  by  meems  of  a servo  con- 
trolled pressure  system.  The  pressure  auid  the  inner  and  outer 
radial  displacements  were  continually  recorded. 

Tests  were  run  at  a constant  temperature  in  the  rubbery 

region,  +26®C,  and  at  two  temperatures  in  the  transition  range, 

-45“C,  -75“C.  The  onset  of  the  glassy  region  for  this  material 

was  determined  to  be  -85 '*C.  Tests  were  conducted  at  the  follow- 

-4  -3  -2 

ing  frequencies:  4 x 10  Hz,  2.5  x 10  Hz,  8 x 10  Hz,  2md  at 
two  pressure  amplitudes.  These  loading  conditions  were  suffi- 
ciently slow  so  the  assunption  that  inertia  effects  are  adssent, 
i.e.  the  deformation  is  quasi-static,  was  effectively  realized 
in  the  experiments. 

The  tensile  specimens  were  cast  as  straight  lengths  of 
synthetic  rubber  100  mm  long  having  square  cross  sections  with 
10  ran  sides.  End  tad>s  were  bonded  directly  to  the  specimens  so 
the  actual  length  could  be  taken  as  the  effective  gauge  length 
and  the  extension  of  the  machine  crosshead  would  determine  the  axial 
extension  ratio  Ai.  Some  measurements  were  also  taken  of  the 
change  in  cross  section  dimensions  which  is  directly  related  to 
Xa.  Tensile  tests  were  conducted  in  an  Znstron  machine  at  the 
sane  temperatures  as  the  sphere  tests  and  at  crosshead  extension 
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rates  of  0.02,  2.0  euid  20  cm/min.  Division  by  the  initial  length 
of  10  cm.  would  give  the  initial  axial  deformation  rates. 

In  110],  it  was  reported  that  two  terms  in  the  series 
expression  for  Eq.  (13),  were  necessary  to  obtain  good 

representation  for  the  two  decades  of  frequency  of  the  sphere 
tests.  This  led  to  five  material  const£uits  in  the  anelastic 
stress  term,  Eq.  (14),  neunely,  ai , (Xz , gi,  gz  and  f at  each 
ten^erature.  In  addition,  the  usual  two  elastic  constants  have 
to  be  specified.  These  seven  constants  were  chosen  to  give  good 
correspondence  between  the  experimental  amd  calculated  results . 
The  same  material  constants  also  led  to  accurate  predictions  of 
uniaxial  stress-strain  curves  with  the  modification  that  the 
parauneter  62  in  the  expression  had  to  be  non-zero  to 
match  the  large  strain  response.  As  mentioned  previously,  a 
non-zero  32  did  not  influence  the  calculations  for  the  sphere 
problem  which  were  limited  to  relatively  low  strains  (ad)out  5%)  . 
In  the  rubbery  region,  +26 “C,  it  was  found  that  a non-zero 
coupling  factor  y,  Eqs.  (8)  emd  (9),  was  more  suitedsle  to  obtain 
good  correspondence.  As  mentioned  previously,  the  shear  modulus 
Uo  was  ta)cen  to  have  a constant  value  at  all  the  test  tempera- 
tures while  the  Leun^  coefficient  X was  taken  to  be  the  same  at 
26*C  and  -45*C  but  differed  at  -75®C  near  the  glass  transition 
temperature.  The  values  of  the  various  material  constants  that 
lead  to  good  correspondence  between  the  test  and  calculated 
results  for  both  the  sphere  and  uniaxial  geometries  are  listed 
in  Table  1. 

Results  for  the  pressure-inner  radial  displacement  relation 
for  the  thick  walled  sphere  at  -75®C  are  shown  in  Fig.  1 for  the 
throe  frequencies.  The  experimental  results  are  shown  by  the 
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continuous  curves  while  the  calculated  points  are  Indicated  by 

symbols.  Other  results  for  the  sphere  problem  are  given  in 

Figs.  2,3,4  at  the  other  test  temperatures.  The  tests  at  +26 “C 

in  the  rubbery  region.  Fig.  4,  show  the  essential  etbsence  of 

anelastic  effects  at  that  temperature.  Unusual  behavior  was 

-4 

observed  at  -45®C  at  a frequency  of  4 x 10  Hz,  Fig.  3,  which 
was  corroborated  by  the  calculations. 

Experimental  uniaxial  stress-strain  results  at  -75 °C  emd 
-45“C  are  given  in  Figs.  5 emd  6.  On  these  graphs,  the  stress 
is  the  "true"  stress  based  on  the  calculated  current  cross- 
sectional  area  and  the  indicated  strain  is  the  Cauchy 
axial  strain  or  Ai  -1.  The  loading  part  of  the  curves  are  shown 
in  Figs.  5 and  6 while  Fig.  7 shows  a more  detailed  picture  of 
the  calculated  loading  and  \inloading  paths  at  -75 "C  for  strains 
up  to  10%.  Calculations  were  also  performed  at  the  -75 "C  con- 
dition for  the  change  of  Poisson's  Ratio  v with  strain  at  two 
loading  rates.  Fig.  8,  where  v is  defined  as 


It  is  noted  that  v decreases  with  both  increasing  deformation 
rate  and  strain.  The  value  of  v at  zero^oa^hg  rate  and  strain 
would  be  given  by  the  standard  elastic  r^l^lon  v « A/2(X+uo) 
which,  for  the  constants  listed  in  Table  1 at  -75"C,  would  give 
V « 0.362. 
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Discussion 

The  computed  deformation  response  for  the  spherical  euid  uni- 
axial geometries  based  on  the  analytical  formulation  showed 
fairly  good  agreement  with  the  experimental  results  over  the 
reuige  of  frequencies  and  temperatures  of  the  program  as  shown  in 
Figs.  1-6.  Although  only  some  typical  results  are  shown  in  these 
diagrams,  the  sphere  tests  were  performed  for  two  different  sets 
of  specimen  dimensions  at  two  maximum  pressure  levels  for  each 
set  with  equally  good  correspondence  [18].  It  was  found  that  the 
constitutive  equations  are  capcible  of  reproducing  the  various 
phenomena  that  were  introduced,  namely  nonlinear  rate  dependent 
response  and  energy  dissipation  (hysteresis)  upon  unloading. 

Figs.  1,  2,  and  5 show  rate  dependent  response  for  both  geo- 
metries while  anelasticity  is  indicated  in  the  loading-unloading 
curves  of  Figs.  1,  2,  3,  and  7. 

The  main  type  of  nonlinearity  inherent  in  this  study  is  of 
the  "magnitude"  type  as  defined  in  [1]  and  [2].  Linear  visco- 
elastic behavior  leads  to  the  relation 

a (a^)*  (16) 

where  a is  the  maximum  amplitude  of  the  anelastic  stress  in 
cyclic  loading.  For  the  anelastic  stress  equation,  Eq.  (14), 
zmd  the  material  constants  given  in  Table  1 at  -75"C,  the  value 
calculated  for  the  power  in  Eq.  (16)  for  the  sphere  tests  shown  in 
Fig.  1 would  be  2.3  which  Indicates  magnitude  nonlinearity. 

Another  result  related  to  this  type  of  nonlinearity  is 
obtained  by  coiq}aring  results  or  calculations  based  on  classical 
linear  elasticity  with  the  nonlinear  elastic  formulation  of  the 
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sphere  problem.  At  circumferential  strains  up  to  cUsout  5%  the 
deviation  between  the  linear  cind  nonlinear  results  is  small 
while  at  10%  strain  the  deviation  rises  to  about  20%. 

A further  aspect  of  nonlinear  behavior  is  found  from  the 
time  dependence  of  the  constant  strain  rate  uniaxial  tests.  It 
is  found  (from  the  Isochronic  data  of  Ted>le  2)  that  the  stress 
at  given  times  does  not  depend  linearly  on  the  velocity.  The 
data  in  Table  2 for  tests  at  -75 ®C  indicates  that  the  response 
is  close  to  linear  up  to  0.2  minutes,  after  which  it  becomes 
nonlinear  from  0.4  to  2 minutes  and  subsequently  reapproaches 
linearity  after  ed^out  four  minutes. 

Rate  dependent  response  is  shown  in  both  the  pressurized 
sphere  aund  uniaxial  tension  tests.  In  both  types  of  experiments 
the  effective  modulus  increases  with  loading  frequency  or  axial 
velocity.  The  material  parameters  in  the  constitutive  equa- 
tion responsible  for  this  behavior  are:  ai,  olz  , gt,  gz . These 
were  determined  to  fit  the  pressurized  sphere  tests  and  were 
then  applied  to  the  uniaxial  tension  tests.  Good  matching  is 
achieved  in  both  types  of  loading,  e.g.  Figs.  1,  2,  5 and  6. 

The  damping  loss  factor  is  defined  by  Eq.  (11).  It  is 
interesting  to  note  that  it  does  show  some  frequency  dependence 
due  to  nonlinear  interactions.  The  confuted  values  of  the  loss 
factor  for  the  pressurized  spheres  at  -7S*C  and  at  a peak  pres- 
sure of  10  N/cm*  are:  0.21,  0.28,  and  0.26  at  the  loading  fre- 
quencies of  4 X 10”“,  2.5  X 10"*,  and  8 x 10"*  Hz  respectively. 

These  values  agree  well  with  damping  measurements  of  low  anqpll- 
tude  free  vibration  torsion  tests  carried  out  at  conqparable  fre-  ; 

quencles.  At  a peak  pressure  of  40  N/cm*  the  computed  loss  j 

factors  at  the  same  three  frequencies  are:  0.072,  0.096  and  0.088  I 

\ 
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Good  matching  of  the  uniaxial  experimental  results  up  to 
eJsout  40%  strain  at  -75®C  euid  100%  at  -45“C  could  be  obtained 
through  the  modification  of  the  elastic  distortion  energy  W 

s 

by  including  a non-zero  value  for  the  62  coefficient.  With 
82=0,  matching  is  limited  to  strains  of  only  15%  at  -75®C.  It 
seems  that  the  restriction  to  fitting  results  with  the  more 
general  expression  for  W^,  Eqs.  (6)  and  (7) , lies  in  its  mono- 
tonic nature.  The  material  stress-strain  curve  at  -75“C  is 
characterized  by  £ui  inflexion  point  which  could  be  accomodated 
if  the  elastic  shear  energy  function,  W , consisted  of  poly- 
nomials  higher  them  the  second.  It  could  then  allow  matching 
of  stress-strain  curves  having  inflexion  points  as  well.  Such 
functions  will,  however,  require  a large  number  of  coefficients 
which  will  create  con^utational  emd  experimental  problems  in 
their  adequate  determination.  It  therefore  seems  that  inflexion 
points  of  the  stress-strain  curve  in^>ose  the  practical  limit 
of  application  of  the  present  theoretical  model. 

The  uncoupled  strain  energy  form  was  found  to  result  in 
good  matching  of  experimental  emd  calculated  results  in  the 
tramsltlon  region  (-45*C,  -75*C) . However,  in  order  to  obtain 
good  matching  of  the  detailed  test  results  for  both  sets  of 
sphere  specimens  emd  for  the  uniaxial  specimens  in  the  rubbery 
region,  26®C,  the  coupled  energy  form  had  to  be  used.  At  that 
temperature,  the  pressure-radial  displacement  curves  of  the 
thinner  walled  spheres  (75mm  - 100mm  diameter)  and  the  uniaxial 
stress-strain  curves  were  extremely  steep  at  the  origin.  This 
behavior  can  be  predicted  from  the  calculations  by  the  use  of 
the  coupled  form  of  the  strain  energy  function. 
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In  conclusion,  the  overall  results  show  that  the  representa-  ] 

I 

tion  capedaility  of  the  analytical  model  is  fairly  good  for  most  j 

engineering  purposes.  The  mathematical  formulation  is  relatively  | 

I 

simple  and  ened^les  accurate  determination  of  the  material  para- 

• I 

I meters.  Since  the  constitutive  equations  are  in  differential  , 

i form,  they  can  be  readily  incorporated  into  numerical  procedures 

I ^ 

• for  solving  stress  analysis  problems  by  step  by  step  methods. 

I The  results  obtained  for  both  the  spherical  emd  uniaxial  geo- 
metries encourage  the  continuation  of  testing  and  evaluation  of  , 

the  theoretical  model. 
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